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Abstract 

Utilizing a complete Lorentz-covariant and local-gauge-invariant formulation, we discuss graphene response 
to arbitrary external electric field. The relation, which is called as Kramers-Kronig relation in the paper, 
between imaginary part and real part of ac conductivity is given. We point out there exists an ambiguity in 
the conductivity computing, attributed to the wick behavior at ultraviolet vicinity. We argue that to study 
electrical response of graphene completely, non-perturbational contribution should be considered. 
PACS numbers: 71.10.-w,72.10.Bg,73.63.-b 



1 Introduction 

Graphene, a flat monolayer of carbon atoms tightly 
packed into a two-dimensional honeycomb lattice, has 
spawned many theoretical and experimental focuses. As 
stated in ref. [1] , graphene plays the role of bridge be- 
tween condensed physics and high energy physics. This 
is attributed to massless Dirac fermion behavior of quasi- 
electron in graphene, i.e. we can treat the quasielectrons 
in graphene as ultimately relativity particles. 

Such behavior also arouses many unusual properties 
of graphene, such as ac and dc conductivity. Many 
attentions '^"^1 are on such topic. However, there ex- 
ist discrepancies in the problem, including the discrep- 
ancy between theories and experiments on dc conduc- 
tivity, the famous missing "i factor", and confliction 
among different theoretical calculations . So far dif- 
ferent theories are almost based on the perturbational 
approximation, even calculations may be did by multi- 
loop diagrams . 

We introduce a correlation function with respect to 
only one variable, the invariant amplitude of spatial-time 
position, X, to study the graphene conductivity non- 
perturbationally. From the correlation function, we show 
there is relation between imaginary part and real part of 
ac conductivity. The function is very close to spectral 
function and we find that the perturbational calcula- 
tions to conductivity only include contributions from free 
valence-conduction electron pairs. Therefore, besides 
these contributions, to compute conductivity completely 
we should also consider other ones, such as excitation or 
impurity. To check the statement, we perform a pertur- 
bational calculation of dc conductivity using the quan- 
tum field theory. This technique guarantees that the for- 



mulation is Lorentz-covariant and local-gauge- invariant. 
We point out that there exists discrepancy among differ- 
ent theoretical calculations, attributed to bad behaviors 
at ultraviolet vicinity of (5-functions. 

We organize the paper as following: In section 2, 
we discuss the electrical response to arbitrary external 
electric field. A discussion on obtaining dc conductivity 
utilizing Kubo theory is also given. The relation be- 
tween imaginary part and real part of ac conductivity is 
listed in section 3. We show an explicit perturbational 
computing for conductivity in section 4 and a brief dis- 
cussion in section 5. 

2 The conductivity under arbi- 
trary external fields 

To perform the calculation wc first give the second 
quantization on graphene briefly. Lagrangian density is 



C = ?A(i7^(9^ - m)i/'. 



(1) 



^ = ^^7o, m is the mass of quasiparti- 



where d^^ — 

cle (To clarify we here endow the quasielectrons with a 
nonzero mass), and 7's are 

70 = = T3, 71 = Pti, 72 = /3t2, 

where ti, T2 and T3 are three Pauli matrices. In the pa- 
per, the repeated indices are generally summed, unless 
otherwise indicated. Furthermore, h = vp = e = I are 
always set. 

The Hamiltonian is then 



(2) 
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Here, Latin indices i, j, etc. generally run over two spa- 
tial coordinate labels, 1,2. (While Greek indices fj,, and 
so on, three spatial-time coordinate labels 0, 1, 2 with 
ccq the time coordinate.) 

Under the second quantization, we have 



where, just as pointed out above, J^{x) is in norm form: 



(27r)2V2^ 



[apu{p)e'^P'^ +blv{p)e^P% 



where po — \/to^"+p^ > q. Solutions to positive energy 
u{p) and to negative energy v{p) satisfy respectively 

v)u — 2po, uu — 2m, uu = ]^ + m, 

v^'v — 2po, vv = —2m, vv — ]^ ~ m. (3) 

The explicit forms of u and v are irrelevant. 
For simplification, we call operator a^{ap), which 
creates(annihilates) quasielectron in conduction 
band, as creation(annihilation) operator which cre- 
ates(annihilates) electron, at the same time, we call 
operator 6j,(6p), which annihilates(creates) electron in 
valance band, as creation(annihilation) operator which 
creates(annihilates) hole. Both electron and hole have 
positive energy, pq. 

Since in the perturbational ground state, the valance 
band is completely filled while the conduction band 
is empty, the en ergy of t he ground state is nonzero, 
Egnd = — J d^P\/m^ + p^ ■ To obtain a Lorentz invari- 
ant ground state, we perform a substraction for all the 
states, E ^ E — Eg„d- Under such subtraction, each 
physics quantity, such as energy, current, etc. should be 
in normal form . 

With the substitution id^^ — > p^, we read Hamilto- 
nian operator eventually 



When concerning of electromagnetic interactions we 
should make a substitution of p^ ^ p^ — eA^ in Eq. 

Denoting = g^^^A'^ with metric matrix g = 
diag{l, — I, —1}, the interacting Lagrangian density is 
^int = —e-ip'yi_iiljA'^ = —J^A'^ and the corresponding 
Hamiltonian is 



^- J d^xC^nt = j d^xJi,{x)At'{x), 



(5) 



d^pd^p' 



= {al,a'pe 



u{p)l^,u{p') 



(271)4^2^ 

+alhl,e^^P+^'>u{p)-f^v{p') 

+bpap,e~'^P+P")='v{p)j^u{p') 

-bl,bpe-''^P-P>v{pMp')}- 

Unlike some papers, we here introduce a factor 
associated with momentum integration, which is at- 
tributed to the Lorentz covariant . 

Both the electron number and the hole number are 
conservative without interaction. However, the only con- 
servation quantity is their difference when interactions 
are included, 



/2» 



/[ d^p 
d^x -.iIj^ := j j;—^{alap-blhp). (6) 



N = 



Generally, the interacting Hamiltonian of graphene 
in external field A^^{xo,x) is described by Eq. ([5]). The 
density operator is p — po + Sp, where pq is the equilib- 
rium density operator and Sp is the leading order correc- 
tion with respect to the external field. 

In Heisenberg picture we have 



iSp{xo) = [Hint,Po]- 



(7) 



Therefore I^l 
Sp{xo) = 



dx'odx'[Jf,{x'ox'),po]A^'{x'ox'). (8) 



Due to the spatial and time translation invariant, at 
zero temperature, on fixed time a;Q, the current density 
at arbitrary position is 

< J,,(x°) >=Tr(<5p(xSJ)J^(xSj0)) 



F = y d'^ppo{al,ap + bl,bp). (4) where 



dxod'^xT^,^{xox)A''{xo + x^x), (9) 



Tf,^{xox) =i < |[J^(xoa;), J;.(00)]| > 



(10) 



with the notation of the ground state | >. 

Noticing in Eq. ^ and (fTU|) . variable xq is defined 
on (— oo, 0), we expand the range of xq onto (— oo, oo) as 

T^,{x) = te{xn) < \ [j^{o),Mx)] \ > 

+ i0{~x„) < \[Mx),J^{0)]\ >, 



(11) 

where 0{xo) is step function: 0{xo) = 1 for xq > 
0,9{xo) = for Xo < 0. Tensor T^y(x) is vanishing for 
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space- like x, i.e. the support of tensor is time-like 
three-dimensional vector x. 

Our expansion is different to the one in Ref. [2] , where 
T^i/ has only forward term or backward term. 

The conductivity should be independent on the gauge 
transformation. This means that, under a local gauge 
transformation A'^ — A" — d'^ f, where / is ar- 

bitrary function with f{xo = Xq) = f{xo = — oo) = 0, 
the current density in Eq. ^ should be invariant. Inte- 
grating Eq. ([9]) by part we find that this requirement is 
satisfied provided d'^T^^ — everywhere. The statement 
can be proven by the facts: l)charge conversation, i.e. 
9" J, = 0; 2)^^e{xo) = -^di-xo) = So^5{xo); 3) the 
equal time commutation relation [Jo{tx), J^(ty)] = 0. 

Tf^u{x) is written as 

T^,u{x) = - g^,n)n(x), (12) 

where □ = d'^dn. — — ~ Hfa;) is scalar 
function with respect to only one variable, invariant am- 
plitude of three-dimensional spatial-time vector x. After 
defining the Fourier transformation of function f{x) as 
f{q) = j dxdtf{x)e'^'^^ with qx = g^x^, we have, then, 

Tfj^Aq) = {qt,qy - q^gf,i,)Tl{q^), (13) 

where n(<7^) is the only function with respect to invariant 
amplitude of q. 

From Eq. ^ it seems that is time dependent 
in the time-invariant external electric field. But this is 
not true. It is enough to illustrate it by a special gauge, 
A" = {0,Ex°,0) otA^ = (0,-i;a;o,0), where A is three- 
dimensional potential and E is the external electric field, 
E„ — {E,0). It is easily to see that 

Mx°o) = -E f ^^e-'^°-°5Mi9o'n(go)(.To+x°), 

(14) 

where ^{qo) = n((7Q,q = 0). Jq (charge density) and J2 
are both vanishing and only Ji ^ 0: 

Ji{xl) = E j ^9on(<Zo') P dxoXoe^«°-°+^-''e^«°-o 

J 27r 9g -I- iq^e 

^ Ea[xl\ (15) 

where the additional factor e*^^" (e is a positive infinitesi- 
mal) is to guarantee that the external field is introduced 
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adiabatically. When e — > 0, -^tt^. — 1, we can re- 

2 

place ^2_^°^^g by unitary. Furthermore, since ^{q^) is 
even function of go, / = 0. We finally 



have a time-independent conductivity 

a = J ^n(g2) ^ U{xo = 0, q - 0), (16) 

where, to obtain a meaningful quantity, we should per- 
form a substraction of 11, i.e. we make a substitution: 
n(2;o = 0, q = 0) ^ n{xQ = 0, q = 0)-n(a;o = -00, q = 
0). In Fourier space, this substraction is the substitu- 
tion n((7o,q) n((7o,q) — 11(50 = 0,q). in the paper, 
we always make such substraction for all the physical 
quantities. 

As expected, we obtain a time-independent current 
density for a steady external field. 

It is not difficult to deduce the response to arbitrary 
external fields. Supposing the external electric field is ac 
with frequency u), A — (0, Eq e^'^^° ,0)e'^^° and substitut- 
ing the potential into expression ([9]), we find that only 
x-component of current density is nonvanishing, 

Mxo) = E,e-o f ^n(gg,0) (17) 
J 27r qo- uj + ie 

Since this potential stands for external electric field 
{Ei,E2) — {iu Eoe^'^'-^° , 0) , the complex conductivity is 

I '^^(llO),;^,^^ (18a) 

= l/^^(,o^o)go + /t^^(9o^o)^^•(l8b) 

Obtaining the dc conductivity from Eq. (jl8ap and 
(jlSbp corresponds to the results obtained from famous 
Kubo theory. However, it is not obvious whether we 
can obtain dc conductivity pop from the limit of Eq. 
p8a|l : Firstly, to obtain Eq. (jlSbp from (fT8a|) we need 
not only convergence of all the integrations, such as 

, 2 

/^n(gg,0) ^^^^^°^^.^^ , etc. but also proper subtrac- 
tion of physics quantities. Secondly, when one chooses 
UJ — directly in Eq. (jl8b|l . he will meet an uncom- 
fortable situation: the first term in Eq. (jlSbp is an am- 
biguous ^. To obtain right dc conductivity we should 
perform computation as follows: we calculate the ac con- 
ductivity in the course of nature from Eq. (fT8|) at ^ 0, 
with proper subtraction. At last we read the dc conduc- 
tivity utilizing the limit of ^ 0. This is just pointed 
out by Kubo which implies that, compare to results 
in references [3], results in references [4] is just the right 
results. Of course, all the results in references [3, 4] are 
obtained by perturbational approach. 

After suitable subtraction, the ac conductivity is 

a=/^n(gg,0) (19) 
J 2tt qo - uj + le 

Generally, n(gQ,0) is not convergent or well defined, as 
will be shown by perturbative calculation in section 4. 
This is relevant to the wick definition of 5— function. 
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Therefore, in order to obtain meaningful physical re- 
sult, we need to perforin subtraction to cancel diver- 
gence. For instance, in Ref. [6] the author has proposed 
a soft (5— function. The subtraction should meet some 
physical criteria. For instance, as found in above para- 
graph, to get the result in Eq. from Eq. p^ . after 
the subtraction 11 is still the function of rather than 
the function of go in Fourier space. In section 4 we shall 
show a explicit subtraction to n(a;o, q). 



3 Kramers-Kronig relation of 
graphene conductivity 

In this section we show a relation between imaginary 
part and real part of graphene conductivity. 

We first give a non-perturbational proof that 11 (g^) 
is real, ^{q^) is real at < obviously, the only needed 
is to prove that n(g^) is also real at q^ > 0. 

At q^ > 0, after inserting complete intermediate 
states J2 |r >< r|, we have for T = Tl^{q), 



d^pe 



-[0{~xa)e{pa) + e{x^)e{-p^)]5{p' - 1) 



(27r)3 p2 



we finally get 



f 

Jo 



dts{t) 



q^-t 



(22) 



with identity ~ -F*i'^(/)- Since s{t) > is real, 
T{q) and therefore H{q) are both real. The spectral den- 
sity is in fact a very important function, which will be 
studied elsewhere '^^l. 

Thus, from Eq. (jl9[) . we obtain an important relation 
between real part and imaginary part of conductivity a 



— 



ds Imf7(s) 



T = 



d^xd^pe^'^^ 
(2^)' 



i{e{xa)-e{-xa)){e'P--e-'ns{p)e{po). 



where the spectral function s{p) is defined as 

2^s{p) =J2< k^(0)|r X r| J^(0)| > {2TT)'S'{p~pr). 

(20) 

The spectral function s{p), which is very close to 
state density, includes not only perturbational contri- 
butions, but also non-perturbational contributions. To 
study non-perturbational contributions, one should con- 
sider, for instance, excitations. 

Inserting dtS{p'^ — t) = 1 at p'^ > 0, one obtains 

T = ^ dts{t) J d^xe''"'i{0{xo) - e{-xo))I, (21) 



where 



^9{po)S{p'-t){e^P--e-^n 
(27r)2 

^e-^P^5ip'-tm-p„)-e{po)). 



Utilizing 



z / mxo)eipo) + 9i~xo)e{-pomp^ - 1) 



{2TTf p^ - t + ie' 



(23) 

This relationship between real part and imaginary part 
of conductivity is beyond the perturbation approach 
and can be considered as Kramers-Kronig relation of 
graphene conductivity. We hope the advanced study of 
graphene may check this relation. 

Eq. is one of the main results of the paper. It 
points out that the electrical response of graphene can 
never be considered as a pure resistance, but a resis- 
tance parallel connected with a capacitor with capacitiv- 
ity n([jj^, 0)/2. Furthermore, due to the obvious relation 
between Imcr and 11, Imcr reflects the state structure of 
graphene. Imcr is therefore a non-perturbational probe 
to detect the state structure of graphene. In this view- 
point, Imcr is a more basic quantity than Recr. Further- 
more, Eq. (|23p is irrelevant to the idiographic interac- 
tions, which means that, the equation holds under very 
general conditions, such as the existence of impurities or 
excitations in graphene. 

More present works reveal that the graphene is rarely 
flat, i.e. there are always ripples in graphene. The non- 
vanishing curvature, raised by ripples, will lead two main 
effects: altering group velocity of quasiparticlc and in- 
troducing effective gauge fields. The first effect possibly 
makes a global correction to conductivity, which may be 
absorbed into the redefinition of spectral function, 11. 
Furthermore, since the holding of Eq. (j23[) is irrelevant 
to idiographic interactions, we conclude that Kramers- 
Kronig relation is still valid for corrugated graphene. 

Eq. ([23]) supplies one possible way to study the dis- 
crepancy of dc conductivity between theories and experi- 
ments. One may first perform perturbational computing 
to imaginary part of ac conductivity and then compares 
the perturbational result to experiments at different fre- 
quency. The discrepancy between these results reveals 
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contribution which can not be ascribed to perturbational 
theory. Furthermore, the complete contributions to dc 
conductivity are not only from the Dirac nodal point, 
but from the spectral structure of carriers. 



The nonzero contribution to K{qo) is jgol > 2m. To sim- 
plify we let m = 0. Thus T{xQ,p = 0) = jiOhxo) - 
eixo))S'{xo). Since n(go, 9 = 0) = ~^Tiqo, q = 0), we 
find the dc conductivity of graphene 



4 Perturbational Calculation of 
dc conductivity 

We here deduce the dc conductivity in perturbational 
approximation. After that we shall discuss an ambiguity 
besides the one pointed out in Ref. [6]. 

We begin the game by a perturbational calculation 
of T^,. Noticing < 0|6p,ap^at^6|,,|0 {2^Y5^[pf - 

Pi)^'^{p'f ~ Pi) f-iid the normal ordering of operators 
in current density, we have, for perturbational ground 
state, 

<0|J^(x)J^(0)|0>= / ,/'rf^,\ f,^e-(P+^>,(24) 

where — v[p')^^^u[p)u{p)^^v[p'). Taking advantage 
of Eq. ([5]), one finds, 



= -2p • p' - 6m^ = -{p + p'Y - 4m 
Tj^ = T is given by a direct computing 



(25) 



T{x) = i{e{xo) - 9{-xo))ia - Am'^)K'{x), (26) 
where 



K'ix) = 



(fip<Pp' 



dp{e 



) f dp{e 



(2^)22po 



tpx p — tpx\ 



(2^)22po 



From 



dp{e 



) f dPpe' 



ipx 



(27r)22po 
dp(e*P^ - e-*P^) f (fiper'P'' 



(2^)2 

Fourier transformation of K' is 

K'ip) = - j ^5{{p - qf - m^)5{q^ ~ m')sgn{qo). (27) 

Here p does not need to be on mass shell, i.e. pq — 
y^p^ + rn^ is not needed, if the integrating factor is 
(Pp. We focus on the case of p = O(0r, p is a time- 
like vector). Letting K{x) = (□ — 4m'^)K'{x) and 
K{qo) = if(go,q = 0), we get 



K{qo) 



q^ + 4m' 
4«o 



e{ql-4m^). 



(28) 



dxi 



dx2{e{~X2)'e{x2))5'{x2), (29) 



utilizing Eq. p6|) . Notice that in above equation we 
have made a subtraction n(xo,p) — ^^q^J'^'^^ = at 

Xq —CO. 

However, the functions, such as 6'{x) and 0{x), are 
not well defined. This means that there possibly exists 
ambiguity in Eq. (|29p . This ambiguity is different to the 
one pointed out in Ref. [6]. 

We consider dc conductivity here. First let S{x) be 
the simplest form, Si{x) = for |a;| > | and 6i{x) = 1/a 
for LtI < §. In this case we obtain cti = = f-!- ~ 
utilizing Eq. This is just the result obtained in 

Ref. [4]. Meanwhile, we can also let S(x) be a some- 
what complex forml^l, S2{x) — ^ ■ At this time 

we get (72 



i+7T_j_ ^ Finally, we can also set 



167r 2tt 

Ssix) = ^cosh2(5f^). We find = ^^^^^^^ ^ 
IJ- numeral value of which is in agreement with 

27r 7r 27r ' ^ 

that in ref. [3]. 

To see the physics meaning of Tq in ^3(2;), we write 
out exphcitly: 0{—x) = j^^^t/Tg ■ The role of Tq is some 

like temperature, which means that Tq^ symbolizes the 
disorder. in Si(x) and 77 in S2{x) 1^1 play the similar 
role. Since cti, a2 and (73 are a—, rj— and To-independent 
respectively, we conclude that the dc conductivity is al- 
most temperature-independent near zero temperature, 
although the conductivity value is ambiguous because of 
the wicked behavior of (5-function. This is verified by 
experments '^^l. 

This is a unexpected occasion that the conductiv- 
ity, a physical observable quantity, varies with different 
definitions of (5— function. The ambiguity is associated 
by the different definitions of (5— function at ultraviolet 
region. One may argue that we can eliminate the am- 
biguity by a standard renormalization schedule in quan- 
tum field theory [^-lol, however, this elimination is still 
contributed to the special definition of 5— function at ul- 
traviolet region. We think that the ambiguity implies 
that the dc conductivity of graphene depends on the be- 
havior of quasielectrons at high energy as well as the 
behavior at Dirac nodal point. This is also pointed by 
Kramers-Kronig relation in Eq. Unfortunately, 
linear dispersion relation of quasielectron docs not hold 
at high energy, which means that, different numeral val- 
ues based on linear dispersion and perturbational ap- 
proaches, need corrections. On the other hand, when 
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we study electrical response of graphene, we always per- 
form calculations utilizing diagrams composed by dif- 
ferent Green functions. To include higher corrections, 
we should use loop diagrams. However, since coupling 
g — 2T:e'^ / ehvp is not small, g ^ 1, comparing to leading 
order, the loop corrections can not be ignored. 

One possibly expects that the correction to conduc- 
tivity given above are not large. If this is the case, our 
computations and others I^''*! indicate that about 30% of 
the full conductivity is from the perturbational contri- 
bution. A question is raised, then, where other contribu- 
tions to conductivity come from. A generalized version 
of Eq. (pS)) tells us that, from the definition of state den- 
sity s{q), perturbational contribution to state density is 



27rsf*(g) = 



4m^ 



(30) 



wicked behavior of (5-function, there is ambiguity in 
graphene conductivity calculations. We argue that the 
full perturbational studies need two corrections: one is 
due to the incorrectness of carrier linear dispersion at 
high-energy and the other is higher order correction. Be- 
sides these corrections, however, there is a furthermore 
correction which is nominated as non-perturbational 
corrections in the paper. This correction comes from 
the contribution of excitations, which is attributed to 
electron-electron interactions. 

Authors are very grateful to Dr. M.G. Xia and Dr. 
E.H. Zhang. This work is supported by the Ministry of 
Science and Technology of China through 973 - project 
under grant No. 2002CB613307, the National Natural 
Science Foundation of China under grant No. 50472052 
and No. 60528008. 



at go > 0. 9 function in this equation reveals that, s^* 
only includes the contribution from pairs of free quasi- 
electron and hole. However, since there are complex in- 
teractions between electron and hole, electron and hole 
may be combined into excitations I^"'' , or in other 
words, it is questionable to consider quasielectrons in 
graphene as 2-dimensional electron gas with no inter- 
acting. To study electrical responses completely, one 
must also consider the contribution of excitations(and 
impurities), attributed to Eq. (|^^ . In standard field 
theory it is difficult to study the contribution pertur- 
bationally. We often nominate the contribution as non- 
perturbational one, such as we did in Ref. [13]. Since the 
coupling is large on graphene, such contribution can not 
be ignored when one consider electrical responses. Ap- 
parently, if m is large enough, the nonzero contribution 
from exciton appears before — Am? . We shall discuss 
such contribution elsewhere '^"'^^ . 

5 Discussion 

The relationship between imaginary part, Imcr, and 
real part, Recr, of ac conductivity is given in paper. Imcr 
depends directly on details of state structure and one 
can study state structure from Imcr. We consider it 
as a non-perturbational probe to detect state structure 
of graphene and it is therefore a very important quan- 
tity. Our formulae are Lorenz-covariant and local-gauge- 
invariant. 

We also perform an explicit perturbational calcula- 
tion using quantum field theory. The computing shows 
that the conductivity is mainly manipulated by the 
momentum-energy relation and there is little nexus be- 
tween the conductivity and state density near Dirac 
nodal point. The computing reveals that, due to the 
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